Spectral characteristics of noisy Josephson flux flow osciflator 
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The current- voltage and spectral characteristics of a flux flow oscillator (FFO), based on a long 
Josephson junction, are studied. The investigations are performed in the range of small bias currents 
and magnetic fields, where the FFO radiates a quasi-chaotic signal with extremely large radiation 
linewidth, and the displaced linear slope (DLS) is observed at the current-voltage characteristic. By 
direct numerical simulation of the sine-Gordon equation it is demonstrated that for large lengths of 
the Josephson junction or in the case of finite matching of the FFO with external waveguide system, 
the DLS with extremely large linewidth is transformed into Fiske steps with very narrow linewidth. 
While there is the common belief that the chaotic regime of the FFO is due to excitation of the 
internal oscillation modes in the "soft" fluxon chain, it is demonstrated that this regime is inspired 
by multiple reflections of the traveling waves from junction ends. 

PACS numbers: 



The flux flow oscillator— (FFO) based on a long Nb- 
AlOx-Nb Josephson junction of overlap geometry, gen- 
erating a broadband quasi-chaotic signal, is a good can- 
didate to be a noisy source for variety of applications, 
such as: nonstationary spectroscopy^, calibration of mix- 
ers, such as based on SIS (superconductor-insulator- 
superconductor) Josephson junctions, and calibration of 
superconductive integrated receiver (SIR)'^,*. In the 
presence of an external magnetic field, the Josephson 
fluxons are formed and under the influence of a bias 
current start to move from one end of the junction to 
another one and are radiated at the border of the junc- 
tion. Depending on the value of a magnetic field and a 
bias current supplied to the FFO, three various modes of 
generation can be outlined. At small magnetic fields, of 
order of a critical field at which fluxons start to enter the 
junction, the chaotic generation with very large width of 
a spectral line, which can reach several GHz at frequency 
of generation from 50 to 200 GHz, is observed^. Further, 
with increase in a magnetic field and frequency of gener- 
ation (from 200 to approximately 400 GHz), almost ver- 
tical Fiske steps are seen at the current-voltage charac- 
teristic (IVC)^-^. Here, due to resonant mode, the spec- 
tral linewidth is smaller than 1 MHz, but thus, because 
of the Fiske steps, there is no possibility of smooth fre- 
quency tuning that complicates application of this mode 
for spectroscopic measurements. With further increase 
of the magnetic field and the oscillation frequency in the 
range from 450 to 700 GHz, the continuous flux flow steps 
are observed at IV C^ Here the frequency can contin- 
uously be tuned in the whole range leading to the width 
of a spectral line from 1 to 50 MHz. 

It is widely accepted in the literature that in the range 
of FFO frequency from 50 to 200 GHz the so-called dis- 
placed linear slope (DLS) is observed at current-voltage 
characteristio^)ifiir^, and the generation linewidth in 
this range is extremely broad. It is usually attributed to 
excitation of the internal oscillation modes in the "soft" 
fluxon chain at weak magnetic flelds^,^. However, most of 
experimental and especially theoretical works dealt with 
the case where the matching between FFO and external 



wave-guide system was pure or even absent^ — . 

The aim of the present paper is the investigation of 
an influence of thermal fluctuations, a finite value of RC 
load and a length of the FFO on the spectral line in this 
quasi-chaotic regime of FFO radiation. We demonstrate 
that for sufficiently long junctions or in the case of good 
matching, which can be modeled by RC-load at FFO 
ends, the DLS with broad linewidth is transformed into 
Fiske steps with extremely narrow linewidth. Therefore, 
the appearance of regime of chaotic oscillations can be 
explained by multiple refiections of the traveling waves 
from junction ends. 

It is known that all basic properties of the FFO can be 
described in the frame of the sine-Gordon equation: 



+ a(t>t 



t>xxt + V- sin(0) + rif{x, t) (1) 



where indices t and x denote temporal and spatial deriva- 
tives, (/? is the phase order parameter. Space and time 
are normalized to the Josephson penetration length Aj 
and to the inverse plasma frequency uj~^, respectively, 
a = ujp/ujc is the damping parameter, ujp = ^leld'hC , 
LUc = 2e/c-Rjv/^, Ic is the critical current, C is the JTJ 
capacitance, i?jv is the normal state resistance, /3 is the 
surface loss parameter, 77 is the dc overlap bias current 
density, normalized to the critical current density Jc, and 
?7y (x, t) is the fluctuational current density. If the critical 
current density is fixed and the fluctuations are treated 
as white Gaussian noise with zero mean, its correlation 
function is: {r]f{x,t)rif{x',t')) = 2ajS{x — x')5{t — t'), 
where 7 = It/{ Jc^j) is the dimensionless noise inten- 
sity. It — 2ekT/h is the thermal current, e is the elec- 
tron charge, h is the Planck constant, k is the Boltzmann 
constant and T is the temperature. 

The boundary conditions, that simulate simple RG- 
loads, see Ref.si^"— , have the form: 

0(0, i)x + rLCL0(O, t)^t - CL0(O, t)tt + (2) 
/3riCL0(O,t),tt+/3<^(O,i)xt =r, 
(t){L, t)x -t- rjiCR(l){L, t)xt + CR(f){L, t)tt + (3) 
l3riiCB.4'{L, t)xtt + I3(I){L, t)xt = T. 
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Here T is the normalized magnetic field, and L is the 
dimensionless length of the FFO in units of Aj. The di- 
mensionless capacitances and resistances, cl,r and r^^fl, 
are the FFO RC-load placed at the left (output) and at 
the right (input) ends, respectively. 
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FIG. 1: Current-voltage characteristic of FFO for various 
matching at junction ends for a = 0.1, L = 10, F = 1.7, 
7 = 0.001. Circles - mismatching, vl ~ rn = 100; diamonds 
- good matching at the radiating end, vl = 2, — 100; 
triangles - perfect matching, = = 1. 
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FIG. 2: Spectral density of FFO for various matching of junc- 
tion ends for a = 0.1, L = 10, F = 1.7, 7 = 0.1, 77 = 0.3. 
Circles - mismatching, vl ~ rn = 100; diamonds - good 
matching at the radiating end, tl — 2, ru — 100; triangles - 
perfect matching, = r_R = 1. 

The computer simulations are performed at the same 
parameters, as in work^: a=:0.1, /3=0, the external 
magnetic field F=1.5-1.9, /3 — 0, the junction length 
i = 10 80. At these parameters, and for absence 
of the matching with external waveguide system, e.g., 
tl = tr = 100, cl = cr = 10, for L = 10, the spe- 
cific displaced linear slope (DLS) appears at I-V char- 
acteristics (Fig. [TJ circles) where chaotic generation is 



observed (Fig. [2l circles). It should be noted that, in 
difference with Fig. 5 of Ref.^ we do not observe high- 
frequency noisy oscillations of spectra: we perform en- 
semble averaging of spectra over 500 realizations and the 
only envelope curve remains, which looks very similar 
to the one in Refi^. All curves in Fig. [5] are com- 
puted for rj = 0.3, see Fig. [T] The account of small 
noise in the system does not change the character of the 
chaotic generation: the IVC remains nearly the same in 
the range v = 0.8 1.4, Fig. [U and the spectrum re- 
mains broad enough. However, matching the FFO with 
external waveguide system makes an essential impact on 
the IVC and the spectral density of a generated signal. 
At good enough matching at the radiating end, — 2, 
rR = 100, the series of Fiske steps appear at the IVC 
(Fig. [U diamonds), the width of the spectral line de- 
creases essentially (Fig. [2l diamonds). With further 
matching improvement, — rR = 1 - perfect matching, 
the corresponding IVC merges into the continuous curve 
(which qualitatively looks very similar to the DLS in the 
mismatching case, only the working frequency range is 
larger), see Fig. [1] triangles, and the generation becomes 
quasi-monochromatic (Fig. [5J triangles). In real devices, 
such smoothing of the IVCs can also be due to surface 
losse s"''"'''"''^ and self-pumping efFect^ d^'"'^^ , which, however, 
become significant for larger frequencies only^^. In both 
cases of good and perfect matching the width of the spec- 
tral line is defined by the intensity of thermal noise of the 
Josephson junctior>iiii^. Thus, the quasi-chaotic gener- 
ation predicted and observed in work^ arises due to nu- 
merous refiections of Josephson vortices from junction 
ends. In a real situation where the FFO is matched 
with an external line for effective transfer of radiation, 
the chaotic mode of generation is replaced by the quasi- 
monochromatic one. 
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FIG. 3: I-V characteristics of the FFO at various values of ex- 
ternal magnetic field, and various conditions of the matching; 
L = 10, 7 = 0.1. 

The matching at the junction ends essentially influ- 
ences the working frequency range. This is especially 
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important for smaller magnetic fields, where in the mis- 
matched case the generation step is nearly absent, see 
Fig. Hfor r = 1.5. We note that the DLS v = 0.8^ 1.4 is 
very stable to the effect of fluctuations, while the steps in 
the range v — 3^5 are completely smoothed for 7 — 0.1, 
compare with Fig. [T] for 7 = 0.001, the DLS in Fig. 
[3] remains nearly unaffected by noise. Due to match- 
ing, for small magnetic fields F = 1.5 1.9 the spectra 
become very narrow (Fig. |4l circles, diamonds and tri- 
angles for fields r=1.5,F=1.7 and r=1.9, respectively). 
The frequency spacing between the peaks of the spectral 
density is proportional to the voltage at the point of the 
current-voltage characteristic of the junction. Thus, in- 
dependently of the choice of the external magnetic field 
value, for the well-matched case the chaotic generation 
disappears. Therefore, it can be recommended for prac- 
tical applications to place Chebyshev transformers not 
only at the output end^/, but also at the input end of 
the FFO to improve matching and to consequently sup- 
press reflections of traveling waves from junction ends. 
This can lead to significant smoothing of IV curves and 
to continuous frequency tuning in all working range from 
50 to 750 GHz. 
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FIG. 4: Power spectral density of the FFO at various external 
magnetic fields F = 1.5, F = 1.7 and F = 1.9; L = 10, 7 = 0.1, 
rj = 0.45. 

Comparison of power spectral densities for various 
lengths of the junction L = 10, L = 15 and L = 20, 
in the perfectly-matched case (Fig. [SJ circles, diamonds 
and triangles, respectively) shows qualitative indepen- 
dence from the length of the junction. 

The situation is completely different for the unloaded 
case. As one can see from Fig. [51 F = 1.7, with increase 
of L the chaotic generation for L = 10 transforms to the 
multi-harmonic generation for L = 20, and to the har- 
monic generation with one strongly dominating peak for 
L = 80. Correspondingly, the increase of the FFO length 
also significantly affects the current-voltage characteris- 
tic, increasing the IVC slope, see Fig. [T] 

Thus, the use of the chaotic mode of generation of the 
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FIG. 5: The spectral density of the FFO radiation for 
the perfectly-matched case at both junction ends at various 
lengths of the junction L = 10, L = 15 and L = 20; 7 = 0.1, 
T] = 0.4. 




0.0 0.4 0.8 1.2 1.6 2.0 



CO 

FIG. 6: The spectral density in the mismatching case for 
different lengths of the junction L = 10, L = 20 and L — 80; 
7 = 0.01, V = 0.25. 



FFO with a broad spectrum has a number of restric- 
tions that does not give the opportunity of its use for 
creation of subTHz non-stationary spectrometers. How- 
ever, other modes of the FFO can be used, where quasi- 
monochromatic generation with broad linewidth is also 
observed (for example, in more important range of fre- 
quencies from 450 to 700 GHz), besides for additional 
broadening of the spectrum an additional stochastiza- 
tion of bias current and magnetic field of the FFO can 
be used. 

In the present paper it is shown that the power spec- 
tral density of the FFO essentially depends on the junc- 
tion length and conditions of the RC-load matching at 
the ends of the Josephson junction. The special case of 
extremely broadband chaotic generation of the Joseph- 
son flux flow oscillator is observed only at the mismatch- 
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ing of the FFO with external waveguide system and 
small junction lengths. Therefore, the appearance of 
regime of chaotic oscillations can be explained by multi- 
ple reflections of the traveling waves from junction ends 
rather than simply by excitation of the internal oscilla- 
tion modes in the "soft" fluxon chain at weak magnetic 
fields^^. 



FIG. 7: The I-V characteristic at various lengths of the junc- 
tion L — 10, L = 20, L = 80, mismatching case; 7 — 0.01. 
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